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Abstract. In this paper wc continue the analytical study of the sabra shell 
model of energy turbulent cascade initiated in [11]. We prove the global exis- 
tence of weak solutions of the inviscid sabra shell model, and show that these 
solutions are unique for some short interval of time. In addition, we prove 
that the solutions conserve the energy, provided that the components of the 
solution satisfy \u n \ < Ck n 1/ ' 3 ( v / nlog(n + 1)) — 1 , for some positive absolute 
constant C, which is the analogue of the Onsager's conjecture for the Euler's 
equations. Moreover, we give a Beal-Kato-Majda type criterion for the blow- 
up of solutions of the inviscid sabra shell model and show the global regularity 
of the solutions in the "two-dimensional" parameters regime. 



1. Introduction 

Shell models of turbulence have attracted interest as useful phenomenological 
models that retain certain features of the Navier-Stokes and Euler equations. Their 
central computational advantage is the parameterization of the fluctuation of a 
turbulent field in an octave of wave numbers A n < \k n \ < X n+1 by very few repre- 
sentative variables. This range of wave numbers is called a shell and the variables 
retained are called shell variables. Like in the Fourier representation of Navier- 
Stokes equations (NSE) , the time evolution of the shell variables is governed by an 
infinite system of coupled ordinary differential equations with quadratic nonlinear- 
ities, with forcing applied to the large scales and viscous dissipation effecting the 
smaller ones. Because of the very reduced number of interactions in each octave 
of wave numbers, the shell models are drastic modification of the original NSE in 
Fourier space. 

The main objective of this work is to investigate the question of existence, unique- 
ness and regularity of solutions of the inviscid sabra shell model of turbulence. This 
model was introduced in [26] and its viscous version was studied analytically in [11]. 
It is worth noting that the results of this article apply equally well to the well-known 
Gledzer-Okhitani-Yamada (GOY) shell model, introduced in [31]. For other shell 
models see, e.g., [7], [18], [20], [31]. A recent review of the subject emphasizing the 
applications of the shell models to the study of the energy-cascade mechanism in 
turbulence can be found in [6]. 

The sabra shell model of turbulence describes the evolution of complex Fourier- 
like components of a scalar velocity field denoted by u n . The associated one- 
dimensional wavenumbers are denoted by fc„, where the discrete index n is referred 
to as the "shell index" . The equations of motion of the viscous sabra shell model 
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of turbulence have the following form 

= i{ak n+1 u n+2 u* n+1 + bk n u n+1 u* l _ 1 - cfc„_iU„_i«„_ 2 ) - f fc 2 u„ + /„, (1) 

for n = 1,2,3,..., and the boundary conditions are u_i = Mo = 0. The wave 
numbers k n are taken to be 

k n = fc A", (2) 
with A > 1 being the shell spacing parameter, and fco > 0. Although the equation 
does not capture any geometry, we will consider L = k^ 1 as a fixed typical length 
scale of the model. In an analogy to the NSE v > represents a kinematic viscosity 
and /„ are the Fourier components of the forcing. 

The choice of the non-linear term in the equation of the sabra model (1) which 
contains only the local interaction between the shells, can be justified in the context 
of the Kolmogorov theory of homogeneous turbulence (see [16], [20]). The theory 
states that there is no interchange of energy between the modes of the velocity field 
with wavenumbers separated at least by "an order of magnitude" . 

The three parameters of the model a, b and c are real. In order for the sabra 
shell model to be a system of the hydrodynamic type we require that in the inviscid 
(u = 0) and unforced (/„ = 0, n = 1, 2, 3, . . . ) case the model will have at least one 
quadratic invariant. Requiring conservation of the energy 



K\ 2 (3) 

leads to the following relation between the parameters of the model, which we will 
refer to it as the energy conservation condition 

a + b + c = 0. (4) 

Moreover, in the inviscid and unforced case the model possesses (formally) another 
quadratic invariant 

OO / \ 

a 



(5) 



The "physically" relevant range of parameters is \a/c\ > 1 (see [26] for details). 
For — 1 < | < the quantity W is not sign-definite and therefore it is common to 
associate it with the helicity - in an analogy to the 3D turbulence. In that regime 
we can rewrite the relation (5) in the form 

oo 

W = £(-l)"fc>„| 2 , (6) 

n=l 

for 

a = l°gA - • (7) 

c 

We call parameters regime corresponding to < - < 1 the 2D regime. This is 
because in that case the second conserved quadratic quantity W is non-negative 
and can be identified with the enstrophy in 2D turbulent flows. We can rewrite the 
expression (5) in the form 

oo 

W-^fc^KI 2 , (8) 

n=l 

where a is also defined by the equation (7). 
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For the parameters satisfying - = — A the sabra shell model becomes "purely 
three-dimensional" , where the quantity (5) scales like the helicity in the 3D Navier- 
Stokes turbulence. It was found (see [13]) that in that case the energy spectrum in 
the inertial range of the GOY shell model (this is also true for the sabra model) has 
the traditional Kolmogorov scaling law /c« 5 ^ 3 . Moreover, while the parameters of 
the model satisfy - = A 2 , the quantity (6) scales like the enstrophy in the Navier- 
Stokcs 2D turbulence. Therefore, this parameters values are usually referred as 
the "purely two-dimensional" regime. In that case the energy spectrum of the 
sabra shell model (see [17], [13]) obeys the scaling law k~ 3 , which is exactly the 
Kraichnan's law of the 2D developed turbulence. 

The famous question of global well-posedness of the 3D Navier-Stokcs and Eulcr 
equations is a major open problem. In [11] we showed global regularity of weak 
and strong solutions of (1) and smooth dependence on the initial data for the case 
v > 0. In this work we address the question of existence of regular solutions of the 
inviscid (v = 0) sabra shell model (1). First, we prove the global in time existence 
of weak solutions with finite energy Similar results for the inviscid GOY shell 
model were obtained recently in ([2]). The existence of weak solutions for inviscid 
hydrodynamic equations with only energy conservation is not known. The only 
other analogue of the 3D Eulcr equations known to posses weak solutions of such 
type is the inviscid surface quasi-geostrophic equation (see [9] and [33]). Next, we 
show that every weak solution u(t)(ui(t), u%(t), . . . ) conserve the energy provided 
that the components of the solution satisfy the decay estimate 

\u n \<Ck^ 3 (V^log(n+l))-\ 

for some positive absolute constant C, namely, provided it is regular enough. A 
similar result for the solutions of Euler equations is known as the Onsager's con- 
jecture (see [32]) and it was proved in [10] (see also [14], [15]). We also give the 
criteria for the weak solutions to remain unique in certain regularity class. 

Next, we show that if the initial data is sufficiently smooth, then the weak so- 
lutions are smooth and unique for a short period of time. Similar results were 
obtained in the context of other discrete models of Euler equations (see, e.g., [21], 
[23], [24], [35]). The well known Beale-Kato-Majda theorem (see [5], [28]) gives a 
criterion for the blow-up of the initially smooth solutions of the 3D Euler equa- 
tions. In section 5 we establish a similar criterion for the inviscid sabra shell model 
equations and use it to show the global regularity of the solutions of the model in 
the 2D parameters regime. This picture is consistent with what is known about the 
global regularity of solutions of the 2D Euler equations (see, e.g., [3], [4], [28], [29], 
[36], and references therein). 

Analytic and numerical study of loss of regularity of solutions of the inviscid 
sabra shell model of turbulence, as well as the dissipation anomaly phenomena in 
that model, is the subject of ongoing research ([12]). 

2. Preliminaries and Functional Setting 

We repeat here for the sake of self consistency the functional settings introduced 
in section 2 of [11]. In particular, the Proposition 1 is a slightly more generalized 
version of the Proposition 1 of [11]. 

Following the classical treatment of the NSE and Euler equations, and in order to 
simplify the notation we are going to write the system (1) in the following functional 
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form 

— + vAu + B(u, u) = f (9a) 

u(0) = u m , (9b) 

in a Hilbert space H. The linear operator A as well as the bilinear operator B will 
be defined below. In our case, the space H will be the sequences space £ 2 over the 
field of complex numbers C. For every u,v € H, the scalar product (•, •) and the 
corresponding norm | • | are defined as 

1/2 



(u,v) = ^2u n v*, \u\ = ^^|u„| 



n— 1 x n— 1 

We denote by {<pj}j^\ the standard canonical orthonormal basis of H, i.e. all the 
entries of <j>j are zero except at the place j it is equal to 1. 

The linear operator A : D(A) —tHisa, positive definite, diagonal operator 
defined through its action on the elements of the canonical basis of H by 

A4>j = k^j, 

where the eigenvalues fc| satisfy the equation (2). The space 

oo 

D(A) = {ueH : \Au\ 2 = ]T k^\u n \ 2 < oo}, 

71=1 

is the domain of A and is a dense subset of H. Moreover, it is a Hilbert space, 
when equipped with the graph norm 

\W\d(A) = \M> VueD(A). 

Using the fact that A is a positive definite operator, we can define the powers A s 
of A for every set 

Vu = (ui,tt2,ti3, ...)) A s u = (kf s ui,k 2s u 2 ,kl 8 u 3 , . . .). 

Furthermore, we define the spaces 

oo 

V s :=D{A s / 2 ) = {u={u u u 2 ,u z ,...) : £ kf \u \ 2 < oo}, (10) 

3=1 

which arc Hilbert spaces equipped with the scalar product 

(«, v) a = (A s / 2 u, A s ' 2 v), Vu, v e D(A S / 2 ), 

and the norm |u|^ = (u,u) a , for every u € D(A S / 2 ). Clearly 

V.CV = HC Vs, Vs > 0. 

Note that the dual space of V s , for every s € R, is V^V- S . We denote the action of 
the element in u £ V- s on v £ V s by 

oo 

(«, V ) a = (A-/ 2 «,^ /2 «) = E u "<- 

n=l 

The case of s = 1 is of a special interest for us. We denote VD{A 1 / 2 ) a Hilbert 
space equipped with a scalar product and norm 

((u,v)) = (A 1 / 2 u,A 1 / 2 v), \\u\\ 2 = ((u,v)), 
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for every u, v € V . The action of the element u € V-\ — V on v € V is denoted by 

oo 

(« lU ) = (r 1 / 2 «,4 1 /M^ )1 X- 

n=l 

Before proceeding and defining the bilinear term B, let us introduce the sequence 
analogue of Sobolev functional spaces. 

Definition 1. For 1 < p < oo and m we define sequence spaces 

i/p 



n=l ' 
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/or 1 < p < oo, 

w m >°° :={«=(« 1 ,u 2 ,...) : ||A m / 2 u||oo = sup (CKI) < oo}. 

l<n<oo 

for w e w m >P we define its norm 

\\u\\ w m., P = \\A m/2 u\\ p , 

where \\-\\ p is the usual norm in the l v sequence space. The special case ofp = 2 and 
m > corresponds to the sequence analogue of the classical Sobolev space, which 
we denote by 

Those spaces are Hilbert with respect to the norm defined above and its corresponding 
inner product. 

The above definition immediately implies that h d = V d , for all d. Moreover, 
Vd C w d '°° and the inclusion map is continuous because 

\\u\\ wd ^ = \\A d / 2 u\\ 00 <\\A d / 2 u\\ 2 = \u\ d . 

The bilinear operator B(u,v) will be defined formally in the following way. Let 
u, v <G H be of the form w^^Li u n4'n an d v — X^^Li v n4>n- Then 



71=1 



B(u, y) = -j'y y l ak n+ iv n+2 u* l+1 + bk n v n+1 u* n _ 1 + 



+ aA;„_iM„_iu„_ 2 + 6fc„_i^„_iu„_ 2 J <j> n , (11) 

where here again u a = u_i = v a = v_i = 0. It is easy to see that our definition of 
B(u,v), together with the energy conservation condition (4) imply that 

B(u, u) = -i ^2 I ak n+1 u n+2 u* n+1 + bk n u n+1 u* l _ 1 ck n —iU n —iU n —2 j (fin ) 
n=l ^ ' 

which is consistent with (1). In [11] we showed that indeed our definition of B(u, v) 
makes sense as an element of H, whenever u € H and v € V or u € and u£ ff. 
The next Proposition is a slightly generalized version of Proposition 1 of [11]. 

Proposition 1. 
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(i) For all d, s,6 G M and for all u £ Vd-e+s, v £ Vd-e-s, and w £ w 1 + 2e >°° 



\(A d B(u,v),w) \<C d , 8 ,g 



\W\\ w i+2e,. 



,\A (d - e+s ^ 2 u\\A (d - 



e-s)/2 



V 



(12) 



where 



C d 



8,6 



= ^|a|(A 1_3<<+3e+S + A -(l-3rf+3e+ 5 ) ) + | 6 |( A 2- + A -(l-3d+30-,)^ _ (13) 

(ii) For all d,s,6 £ M. and for all w £ Vd-e+s, v £ Vd-e-s, and u £ w 1+29 '°° 
\(A d B(u,v),w) d _ e+s \ < c^^llulL^^I^-^/^ll^- 6 ^)/ 2 ^, (14) 

where 

c d ,s,e = (\a\(\-( 2d - 2s ^ + X 2d - 2s ) + \b\(X 1+d - 30 - s + x 1 + d + 38 -')^ . (15) 

(iii) For every d,s £ M, the operator B : V 2 d-2s x V\+2s — * ^2d an d B : 
V 1+2s x V 2d - 2s — » V 2 d is bounded and 



A d B(u,v) 



< 



C dyS - dyS \\u\\ w l+2e,oo\A d S v\, 

C'd, s -dA v \\u ] i+^.~\A d - s u\, 



where the constant Cd, s ,e and Cd, s ,e were defined in (13) and (15). 
(iv) For every u £ Vd, v £ V\- 2d , for all del, 



(B(u,v),u) d = - (B(u,u),v) 



l-2<2 > 



and 



Re (B(v,u),u) d = 0. 
Proof. To prove the inequality (i), we write 



(16) 

(17) 
(18) 



(A d B(u,v),w) 



1+26 



{ak n+1 kl d v n+2 ul +1 w* n + bk 2d+1 v n+1 wlu* n _ 1 

nl ^ 

+ afc^/c„„iw>„_i?;„_2 + &A£ d fc n _itu*t; n _iit n _2 

OO 

1=1 

bX^{-'v n+1 k^w^{ + '<-i 

aA -(l-3«I +M+ .) fc l + 2 e< ^- ?+X _ iA .d- e2 - aVn _ 2 
6 A-( 1 - 3 "+ 3 «-«)fci+ M <fct;-'«n-l^: fl 2 + *«n-2 



< 



+ 



+ 



< 



<c d , M || w |L 1+2 ^|^- e + s )/ 2 u ||A( d - e - s )/ 2 w |. 

In the same way we prove the inequality in (ii). 
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In order to prove the statement (iii) we apply (14) to obtain the first inequality 



A d B{u,v) 



sup 



< 



{A a B{u,v),w) 

< SUp Cd, s -d,s||u|| W l + 2n,oa \w\ \A d ~ S v\ < 

\w\ — l 

< Cd )S _d )S ||ti|| t(; l+2.,oo|A d - S u|. 

The second inequality is proved similarly. 

Finally, the statements (iv) follow directly from the definition of the bilinear 
operator B(u,v), the energy conservation condition (4) and the inequality (12). 

□ 

3. Weak solutions of the inviscid shell model 

Let us consider the inviscid sabra shell model problem 

du 

— + B(u,u) = f (19a) 

u(0) = u m . (19b) 

One of the main properties of the sabra shell model of turbulence is the locality 
of the non-linear interaction. This property allows us to prove the global existence 
of weak solutions with the finite energy to the inviscid shell model in the following 
sense. Similar results in the context of the GOY shell model with the stochastic 
forcing were obtained recently in ([2]). In the rest of the section we give the sufficient 
criteria for the weak solutions to conserve the energy, and investigate the question 
of the uniqueness of weak solutions. 

Definition 2. Let < T < oo, then u(t) G L°°([0,T], H) f]C([0,T], H w ) is called 
a weak solution of the system (19) on the interval [0,T] if for every < t < T it 
satisfies 

(u(t), v)+ f (B(u(s), u(s)), v) ds = (u in ,v) + (/, v) , (20) 
Jo 

for every v ^V. 

Observe that if u(t) = (u\(t), U2{t), Us(t) . . .), then u G C([0, T], H w ) is equiva- 
lent to u n {t) e C([0, T], C), for every n = 1,2,2,, ... . 

Theorem 1. Let u ln , f <G H, then for every < T < oo a weak solution 

u(t) e i oo ([0,T],i/)f|C([0,T],i/ u) ), (21) 

in the sense of Definition 2 exists. In addition, 

^eL°°([o,r],^i)- (22) 

dt 

Proof. Let us fix m > 1. Denote by P m - the orthogonal projection in H onto the 
first m coordinates and Q m = I — P m . The Galcrkin approximating system of order 
m for equation (19) is an m-dimensional system of ordinary differential equations 

du m 

— + P m B(u m , u m ) = P rn f (23a) 
u m (0) = P m u m . (23b) 
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First observe that the nonlinear term of equations (23) is quadratic in u m . There- 
fore, by the theory of ordinary differential equations, the system (23) has a unique 
solution on some finite time interval [0,T^). Let us now take the inner product of 
both sides of the equation (23a) with u m and using subsequently inequality (18) of 
Proposition 1 and the Cauchy-Schwartz inequality we get 

lj t \u m \ 2 = (Pmf,u m ) < \P m .f\\u m \ < (24) 

from which we conclude that 

\u m (t)\ < \u m (0)\ + \f\t< \u m \ + \f\t. (25) 

Therefore, u m is finite in the H norm for all t < oo, hence we can extend the 
solution of the problem (23) to all t £ [0, oo). 

Let us fix < T < oo. Then, from the relation (25) we may conclude that 

sup |<(*)|<C, 

0<t<T 

for some constant C > depending only on u m , /, and T. Moreover, writing the 
equation (23a) in the componentwise form 

«:w = <(o)+| t ! ( a uc 2 (Ci) , +^: + i(Ci)'+ 

+ aAv.-iiCiiCa + bk n - lU ™_ lU ™_ 2 y S + /„. (26) 

For < t < T, we get that for every n there exists a constant C n , independent of 
to, such that 

ll^Cllciao.nc) < c n . 

Applying Arzela-Ascoli theorem we conclude that for every n there exists a subse- 
quence (m£)fc>i such that Un k converges uniformly to some u n , as k — > oo. More- 
over, by a diagonalizing procedure we can choose a sequence (mfc)fc>i, independent 
of n such that u™ converges uniformly to u n G C([0,T],C) and we denote 

u(t) = (u 1 (t),u 2 {t),u s (t), ...). 

Using the uniform convergence it is easy to show, passing to the limit in the ex- 
pression (26), that u(t) satisfies the weak form of the sabra shell model equation in 
the form 

(u(t),v n ) + f (B(u(s),u(s)),v n ) ds = (u m ,v n ) + (f, v n ) , (27) 
Jo 

for every v n £ H with the finite number of components different from zero. 

Moreover, we need to show that u(t) £ L°°([0, T], H). The sequence {u m } m >i 
is uniformly bounded in H (see (25)), and hence 

u m is bounded in every L p ([0, T],H), for 1 < p < oo. (28) 

Therefore, we conclude that there exists a subsequence (mk)k>i such that u mk 
converges to w(t) in the weak-* topology of L°°([0, T], H), and by definition it is 
not hard to see that the limiting function is indeed w(t) = u(t). In addition, by the 
inequality (iii) of Proposition 1 , we get 

|B(u,«)|_i = |yr 1/2 B( M ,w)| < CM 2 , 
for C = C -1/2,0,-1/2 ( scc (13)), concluding that B(u,u) £ L°°([0, T], V-\). 
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Finally, let v € V, and v n — P n v, with the finite number of components being 
different from zero, converging strongly to v. Then letting n — > oo in the rela- 
tion (27) we conclude that u(t) satisfies the equation (19) in the weak sense of 
Definition 2. 

□ 

The next Theorem gives a partial answer to the question: under which conditions 
the weak solutions conserve the energy? 

Theorem 2. Let u(t) be a weak solution, whose existence is proved in Theorem 1, 
satisfying 

u(t)e L°°([0,T\,V 1/3 ), (29) 
for some T > 0. Then for every t G [0, T] 

\u(t)\ 2 = \u in \ 2 + (\f,u{s))ds. (30) 
Jo 

Proof. If a weak solution satisfies (29), then according to the inequality (iv) of 
Proposition 1, the equation (19) can be considered as the equation in the space 
V_i/z. Applying the operator A~ 1/>3 to both sides of equation (19), and taking 
an inner product with A 1//3 u(t) in the space H we get, using the identity (iv) of 
Proposition 1, 

~l«(*)| 2 = (/.«(*)). 
from which the statement follows. □ 

As we already mentioned in the introduction, our result for the sabra shell model 
of turbulence is reminiscent of the Onsager's conjecture for the Euler equations (see 
[10], [14], [15], [32]). However, the criterion given by Theorem 2 is not sharp. It is 
easy to give an example of a solution of the inviscid sabra shell model of turbulence, 
which stays merely in H, but still conserves the energy. To see this, consider the 
forcing f(f 1 ,f 2 ,...), where 

= 1,3,6, ... , 



fn 



\ 0, o/w. 



Then solution u(t) = (ui(t), U2(t), . . . ), where u n (t)-, for n = 1,3,6,..., is a 
weak solution of the sabra shell model, corresponding to the zero initial condition. 
The function u(t) is only in H for every t < oo, however it is easy to see that 
it conserves the energy. Clearly, this example is pathological in a sense that all 
non-linear interactions are absent due to the wide gaps between the excited modes, 
however it shows that the result of Theorem 2 is not sharp. Moreover, it is not 
known when the weak solutions of the inviscid sabra shell model dissipate energy. 
These questions will be studied in the forthcoming work [12]. 

The final result of this section gives criterions for the uniqueness of weak solu- 
tions. 

Theorem 3. 

(i) Let u(t),v(t) be two weak solutions, whose existence is proved in Theo- 
rem 1, satisfying 

u{t),v{t)eL\{0,Tlw l >°°), (31) 
for some T > 0, and u(0) = v(0). Then u{t) = v(t), for all t e [0,T]. 
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(ii) If u(t),v(t) are two weak solutions, satisfying 

u(t) e L\[0, T], w 1 ' 00 ) p| L°°([0, T], V 1/3 ), (32) 

and 

t ) (t)eI co ([0 ) Jl ) V 1/3 ) ) (33) 
for some T > 0, with u(0) = v(0). Then u(t) = v(t), for all t G [0, T]. 

Proof. First, let u, v <G i 1 ([0, T], w 1 ' 00 ) be two weak solutions of the inviscid sabra 
shell model with the same initial conditions. Denote w = u — v satisfying 

dw , , 

— + B(u, w) + B(w, v) = 0, 

with w(0) = 0. Using the fact that u, v and w satisfy (22) and (21), we are allowed 
to apply the operator A^ 1 / 2 to both sides of the last equation and then take the 
inner product of both sides with w in H to conclude 

\j t \ A ~ 1/Aw \ 2 < \(A-^ 2 B(u,w),w)\ + \(A-^B(w,v),w)\ = 
= \(A-^ 2 B(u,w),w)\ + \(B(w,A~ 1/2 w),v)\ < 
KC^MUoo + M^A-V'wf, (34) 

where we subsequently used parts (iv), (i), (ii) of Proposition 1 and 

C\ = C-i/2. -1/2.0 + c -i/2,o.o- Applying the Gronwall's inequality to (34) we get 

lA^^wit)] 2 < |A- 1 / 4 w;(0)| 2 e c ' 1 .''o(ll t '( s )ll™ 1 - + ll u ( s )ll™ 1 -) <is , 

for t € [0,T], concluding the proof of the part (i). 

To prove the part (ii) of the theorem, let u{t) be the solution of the inviscid sabra 
shell model satisfying (32). Let v(t), satisfying (33), be another weak solution with 
the same initial data v(0) = u(0). Note that, in particular, both u(t) and v(t) 
conserve the energy, according to Theorem 2. Denote w = u — v satisfying 

^ + B(u,w) + B(w,u) +B(w,w) = 0, (35) 

with w(0) = 0. Just as in the proof of Theorem 2, we can consider (35) as an 
equation in the space V_i/ 3 . Therefore, applying the operator A^ 1 ^ 3 to both sides 
of equation (19), and taking an inner product with J 4 1 / 3 w(t) in the space H we get, 
using the parts (i) and (iv) of Proposition 1, 

< \{B(w,u), w) 1/3 \ < c 2 |Mlu.i>~M 2 , 

for C2 = Co, o,o- It follows that 

Hi)| 2 < | W (0)| 2 e C2 ^ll«( s )IL^- ds , 
for t € [0, T], finishing the proof of the theorem. □ 
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4. The short-time existence and uniqueness of strong solutions 

The uniqueness of the weak solutions, whose existence was proved in Theorem 1 
is not known. In this section we prove that the weak solutions of the inviscid case 
sabra shell model are unique as long as they stay smooth enough, at least on the 
short-time interval [0, T»), where the time T» depends on the parameters of the 
problem (a, b, c, fc and A), as well as the initial data u m and /. Let us consider the 
inviscid sabra shell model problem (19) as an ordinary differential equation (ODE) 
in the Hilbert space Vd, for d > 1. The main theorem of this section shows the 
short time existence and uniqueness of solutions of the equation (19). 

Theorem 4. Let u m e Vd and f £ Vd for some d > 1. 

(i) There exists a time T > 0, such that the inviscid problem (19) has a unique 
solution u(t) satisfying 

u(t)£C 1 ((-T,T),V d ). 

(ii) Moreover, if f £ V2d-1, then 

^ &C ((-T,T),V 2d -i). 

(iii) The unique solution to the inviscid sabra shell model (19) either exists 
globally in time, or there exists a maximal positive time of existence T* > 
such that 

u(i)eC7 1 ([o,T,),y d ), 

and 

limsup|w(i)|d = oo. 

A similar statement can be formulated for the maximal negative time of 
existence. 

In our notation C 1 ((— T, T), Vd) denotes continuously differentiable functions on 
the interval (—T,T) with values in Vd- Our proof is based on the classical Picard 
theorem for ODEs in Banach spaces (see, for example, [19], [28], [34]). 

Proof. (Of Theorem 4) 

Let us write the system (19) in the form 

^=F(u), u(Q)=u m , (36) 

where F(u) = f — B{u,u). Fix d > 1, then according to the part (iii) of Proposi- 
tion 1 the operator B(u,u) maps V d into V^-i Q Vd- Hence, the mapping F(u) 
maps Vd into itself. Moreover, for every u, v £ Vd we have the following estimates 

\F(u) - F(v)\ d - \B{u, u) - B(v, v)\ d < 

\B{u-v,u)\d + \B(v,u-v)\ d < C 2 (\\u\\ + |H|)|u -«!«,, 

where the last inequality follows from relation (16), and C 2 = C'd.-d,o + c<j,-d.o > 
(see (13)). Therefore, we conclude that the mapping F(u) is locally Lipschitz 
continuous, and we are able to apply the Picard theorem, completing the proof. 

Part (iii) follows by the straightforward application of classical theory of ODEs. 
To prove part (ii) we apply inequality (16) to both sides of the equation (19). 

□ 
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Using the part (iii) of Theorem 4 we will be able in the section 5 to derive a 
criterion for the blow-up of the solutions of the inviscid sabra shell model and to 
prove the global (in time) existence of the unique, regular solutions in the particular 
case of the 2-D regime of the inviscid sabra shell model. 

5. A Beale-Kato-Majda type result 

The Beale-Kato-Majda theorem (sec [5], [28]) states, citing the original article, 
"if a solution of the Euler or Navier-Stokes equations is initially smooth and loses 
its regularity at some later time, then the maximum vorticity necessarily grows 
without bound as the critical time approaches." More precisely, if the initially 
smooth solution of the Euler equations cannot be continued beyond the time T* , 
and T* is the first such time, then 

lim / \\u>(-, s)\\L°°ds = oo, 
t^T-- J 

where lo = curlv is the vorticity and v is the velocity field of the Euler equations. 

Our goal in this section is to derive a similar criterion for the loss of regularity 
of the solutions of the inviscid sabra shell model. Note that in our case, the analog 
of the L°° norm of the vorticity would be the i°° norm of the velocity "derivative" , 
namely 

||tt||^i,oo = sup fc„|u„|. 

l<n<oo 

Clearly, if this quantity becomes infinite at some finite moment of time, then all 
higher norms, namely \u\d, for d > 1, become unbounded at the same time. How- 
ever, in the spirit of the Beale-Kato-Majda result for the Euler equations, we show 
that the opposite is also true. In other words, we show that if a regular solution 
u(t) of the inviscid shell model problem loses its regularity for the first time at the 
time T, then 

rt 

\u(s)\\ w i,oods — ► oo, as t — > T~ . 



o 

For simplicity we would like to focus on the inviscid sabra shell model problem 
(19) without forcing 

C ^ + B(u,u) = (37a) 

dt 

u(0) = u in {x). (37b) 

Theorem 5. Let u in e V d , for some d > 1. Let u(t) e ^([0, T»), V d ) be the 
solution of the inviscid shell model equation (37), where T* is its maximal positive 
time of existence. Then, either T» = oo or 

ft 



and hence 



lim / Hu^ll^i.oodroo, 
t^rr Jo 

limsup||'u(£)|| to i,oo = oo. 



Proof. Let us fix d > 1 and consider u(t) - the unique solution to the sabra shell 
model equation (37). According to the part (ii) of Theorem 4, both sides of the 
equation (37a) lie in the space Vid-\- Therefore, we are allowed to apply the 
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operator A d l 2 to both sides of the equation and take the inner product in H with 
A d l 2 u G H. After taking the real part we obtain 

\l\u\l l = Re(A^B(u,u),A^u) = 

oo , 

= -Im^2 [ ak n+1 k 2 ^u n+2 u n+l u* n + bk 2 n d+1 u n+1 u* n u* n _ 1 + 

n=l ^ 

+ ak^k n -iu* n u n -\u n -2 + 6fc^fc„_iu*M„_iu„_ 2 ^ = 

oo 

= E d ImY / k 2 n d+1 u n+1 u* n u* n _ 1 , (38) 

n=2 

where we denote by 

E d = a(\ 2d - \- 2d ) + b{\ 2d -1). (39) 
Applying Cauchy-Schwarz inequality to (38) we get 

j t \u\ 2 <2\E d \\\u\\ w ^\n\l 

which is valid for t € [0,T], and for every T < T*. Using Gronwall's inequality we 
conclude 

\u(t)\ d < KO)Uel £d IJo T ll«MILi,~^. (40) 
The Theorem follows after letting T — > T~ in the last inequality. 

□ 



6. The "two-dimensional" regime 

Recall that for the range of parameters satisfying < c/a < 1, corresponding to 
the "two-dimensional" regime, the sabra shell model possesses two different positive 
quadratic invariants, one is associated with the energy and the second is associated 
with the enstrophy in the analogy with the 2-D Euler equation (see [17]) 

W= \A d °/ 2 u\ = \u\ do , (41) 

where 

d ° = l l0 ^(l)- ( 42 ) 

In this case E da , defined in (39), equals 0. Moreover, if d > 1, the following 
inequality holds 

||u|Li.°o < |u|do, 

for every u G V do . From the relation (42) we conclude that the condition d > 1 
corresponds to the case when parameters of the inviscid sabra shell model satisfy 

< - < \~ 2 . (43) 
a 

It is well known that the Euler equations of the ideal incompressible fluid in 2-D 
possesses global in time, unique, regular solution (see, for example, [3], [4], [22], [28], 
[29], [36]). The same statement is true for the inviscid shell model of turbulence, 
namely. 



Corollary 1. (Global Existence) Let do be defined by the relation (42). 
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(i) Let the parameters a, c, and A of the inviscid sabra shell model (37) satisfy 

- > A . 

a 

Then for u m € Vd a , there exists a weak solution u(t) to the inviscid problem 
(37) satisfying 

W (i)ei°°((-oo,oo),y do ). 

(ii) The weak solution u(t) conserves the enstrophy (4-1), for all t e [0, T], 
provided 

u(t)eL°°([0,T],V {1 _ 2do)/3 ). 

(iii) // the parameters a, c, and A of the inviscid sabra shell model (37) satisfy 
the relation (43), then for u m G V s , s > do, there exists a unique global 
solution u(t) to the inviscid problem (37) satisfying 

u(t) e ^((-oo, oo), Vi). 

The proof of the part (i) of Corollary 1 is essentially the same as the proof of 
Theorem 1. The proof of the part (i) is similar to that of Theorem 2, and the part 
(iii) follows from the criterion, proved in Theorem 5. 

The comprehensive numerical study of the shell model of turbulence in the 2-D 
parameters regime was performed in [17] for the sabra model, and previously in 
[13] for the GOY model. In particular, it showed that parameters setting defined 
by relation (43) corresponds to the enstrophy and energy equipartition across the 
inertial range. Therefore, our rigorous result on existence of the solutions of the 
inviscid shell model (37) globally in time supports these findings. 

It was also found numerically, that for the parameters satisfying | > A~ 2 , the 
shell models exhibit the direct enstrophy cascade in the inertial range and the 
energy distribution becomes close to the Kraichnan's dimensional prediction 

(KI 2 )~^ i(1+logA(a/c)) , 

with small corrections, for n/ <C n <C nd, where n/ is the largest wavenumber of 
the forcing and nd is the Kraichnan's dissipation wavenumber (see [13]). Note that 
for | = A~ 2 this estimate exactly coincides with the well-known prediction k~ 3 
for the energy spectrum of the 2-D developed turbulence (see [25]). Using these 
reasonings we conclude that for 

- a > A" 2 , (44) 

the inertial range of the sabra shell model with non-zero viscosity will scale like 
(|w„| 2 ) <~ k~ 2+s , for some positive 6. Therefore, it is natural to expect that if the 
viscosity tends to zero, or equivalently, the dissipation scale nd tends to infinity, the 
solutions of the inviscid sabra shell model, for parameters satisfying relation (44), 
will blow-up in finite time, for some initial conditions, according to the criterion 
proved in Theorem 5. 

7. Cubic invariant and Hamiltonian structure 

In practical numerical simulations of the sabra shell model one is limited to 
consider a truncated model of N equations, setting u n = 0, for n = N+l, N + 2, .... 
It was shown in [27] that such a finite system in the inviscid and unforced case 
possesses a Hamiltonian structure for a specific value of the parameters. In this 
section we will state, based on our results on the existence of the solutions of 
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the inviscid sabra shell model, that the infinite system of equations also has a 
Hamiltonian structure. 

By rescaling the time and taking into the account the energy conservation as- 
sumption (4), we will assume that 

= 1, b = — e, c = e — 1. 

Let us fix e = (\/5 — l)/2 - the golden mean satisfying e 2 = 1 — e. In that case we 
can rewrite the equations (1) in the equivalent form 

= ik n +l(u n +2Un+l ~ -^ u n+l u n-l + ^"11-1^-2) - V~k n U n + /„, (45) 

for n = 1,2,3, ... . In that case the inviscid sabra shell model without forcing has, 
fomrally, a cubic invariant of the form 



y^efcpf - - j {u* n+1 u n u n -i+ c.c.), 

n=l V £ / 



(46) 



where c.c. stands for complex conjugate. 

Following the method of [27] , we perform the following change of variables 



and 



an = ~^/2> for even n, 



a n = — —77, for odd n. 



The modified equations then take the form of 

^ = -ik e(\^/e) n ^\yfea n+2 a n+ i + a* +1 a„_i + ^=a* _ 1 a„_ 2 ^ , (47) 
for n = 1, 2, 3, .... Finally, the Hamiltonian takes the form 

oc 

H = Y / ^n, (48) 

71=1 

where 

U n = ek a (Xy/e) n (a n+1 a n a* n _ 1 + c.c). (49) 
In order to see that H is indeed a Hamiltonian we note that 

and the equations of motion (47) satisfy 

da n _ .dU da* n _ . dH 
dt 9a* ' dt da n 

Both the cubic invariant I and the Hamiltonian H are defined by infinite sums. 
Therefore, it is natural to ask when those definitions make sense, namely when the 
sums converge. 

Lemma 1. For A 2 > e _1 the Hamiltonian (4-8) is well defined for all u G V. 
Proof. The Lemma follows by a simple application of Holder inequality. □ 
Finally, we can conclude the following. 
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Corollary 2. Let the parameters of the inviscid sabra shell model (37) satisfy 
a = 1, o = , c= , A > 



VE-i' 



Then the inviscid sabra shell model (37) with initial data in for d > 1, is a 
Hamiltonian system with the Hamiltonian defined by relation (48), as long as a 
solution of the model exists. 



8. Conclusions 

In this work we continued the analytic study of the shell models of turbulence, 
initiated in [11]. We established the global existence of weak solutions and showed 
that strong solutions remain regular and unique for some short period of time. 
Moreover, we showed that the solutions for the "two-dimensional" range of param- 
eters remain regular and unique globally in time. In addition, we established a 
Beale-Kato-Majda type criterion for the blow-up of the initially smooth solutions. 

We showed that for some parameter regime the sabra shell model is an infinite 
dimensional Hamiltonian system. In contrast to the Euler equations, which possess 
a quadratic Hamiltonian function (see, for example, [1], [30]), the Hamiltonian of 
the inviscid sabra shell model is cubic. 

We showed that the weak solution u(t) = (ui(t) 7 u 2 {t) 7 . . . ) conserve the energy 
provided that the components of the solution satisfy 

K\ <Ck^ 3 (V^log(n+l))-\ 

for some positive absolute constant C . A similar result for the Euler equations is 
known as the Onsager's conjecture (see [32]) and it was proved in [10] (see also [14], 
[15]). The question of whether less regular solutions dissipate energy remains open. 

The question of the possible loss of regularity and uniqueness for the initially 
smooth solutions of the inviscid sabra shell model outside of the "two-dimensional" 
range of parameters still remains open. 

The dimensional argument for the viscous [y > 0) GOY shell model, which 
are also applicable to the sabra model, indicate that in the "three dimensional" 
parameters regime — 1 < ^ < the velocity field scales like 

(| Un |)-fcn l(1+l0SA|a/C|) , 

at the inertial range (see [13]). Therefore, at least for the parameters regime satis- 
fying 



> A" 2 , (50) 



we might expect the blow-up of the inviscid sabra (as well as GOY) shell model of 
turbulence. 

We would like to mention that the techniques, used to prove the blow-up for 
other discrete models of Euler equations (see, e.g., [8], [21], [23], [24], [35]) could 
not be applied directly in the case of the sabra shell model of turbulence. The study 
of possible loss of regularity of the inviscid sabra shell model in different parameters 
regime is the subject of ongoing work ([12]). 
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